Backscattering Differential Ghost Imaging in Turbid Media 
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In this Letter we present experimental results concerning the retrieval of images of absorbing 
objects immersed in turbid media via differential ghost imaging (DGI) in a backscattering configu- 
ration. The method has been applied, for the first time to our knowledge, to the imaging of small 
thin black objects located at different depths inside a turbid solution of polystyrene nanospheres and 
its performances assessed via comparison with standard imaging techniques. A simple theoretical 
model capable of describing the basic optics of DGI in turbid media is proposed. 



Ghost imaging (GI) is an optical technique for the re- 
trieval of images via intensity correlation of two corre- 
lated light beams. The first experimental approach and 
theoretical explanation of GI was quantum-like [1], but 
after a long-standing debate [2], it was finally demon- 
strated that GI can be also realized with classical light 
beams [3|,|4| . Thermal GI, for instance, is performed with 
two spatially correlated speckle beams obtained by using 
a rotating ground glass and a beam splitter. The object 
beam illuminates the object and is collected by a bucket 
detector with no spatial resolution, while the reference 
beam is recorded by a spatial-resolving detector, for ex- 
ample by a charge coupled device (CCD) camera. Recent 
improvements of the GI protocol have been achieved via 
computational GI that uses computer controlled spatial 
light modulators [5j, compressive sensing GI where the 
algorithm for the data analysis benefits from the sparsity 
properties of the object [6] and via Differential Ghost 
Imaging (DGI), which has been shown to perform much 
better than conventional GI when imaging weakly ab- 
sorbing objects 0- 

The potentialities of GI with respect to standard (not 
correlated) imaging resides in its ability of forming im- 
ages without necessity of any pixelated detector placed 
nearby the object. Thus GI is a good candidate for imag- 
ing objects immersed in optically harsh or noisy environ- 
ments such as, for example, in a turbid medium or in 
the presence of optical aberrations. Recent applications 
of GI in this direction include imaging in presence of at- 
mospheric turbulence jsj, fluorescent ghost imaging [9| 
and transmission GI in scattering media [10[. All these 
works have raised the very interesting debate whether GI 
is intrinsically more powerful than standard imaging and 
can be used, for example, as a standoff sensing technique 
which is immune from atmospheric turbulence [11] . 

Following this debate, we propose for the first time in 
this letter the use of DGI for the imaging of absorbing 
objects immersed in a turbid medium, in proximity of 
its surface. We adopt a backscattering configuration of 
the bucket light detection similar to the schemes used 
in biomedical tissue imaging and we are able to provide 
information on the transmittance of the object as a func- 



tion of its depth inside the turbid medium. We also com- 
pare the performances of DGI with standard techniques, 
where the imaging is performed with a lens and a CCD, 
showing that the two methods are fairly equivalent. 

The experimental setup for the DGI configuration is 
sketched in Fig. [U The pseudothermal source, operating 
at A = 0.532/im, produces a rectangular collimated beam 
of deep Fresnel speckles [12| with a constant transver- 
sal size S x ~ 82/im and longitudinal coherence length 
5 Z ~ tt^/A ~ 40mm. The beam area is = 46.6mm 2 
and contains iV spe ckie — 6900 speckles of coherence area 
A co h = £ 2 [13[. The reference beam intensity 12(^2 ) 
is recorded at a distance Z2 — 250mm from the source 
by a CCD camera with pixel size 6.67 jam <C S x . The 
intensity ii(xi) hits the object at a distance z\ = Z2 
and is collected with a bucket detection in backscatter- 
ing. The object, characterized by a spatial transmittance 
T(x) over the same area of the beam, is immersed in a 
turbid solution contained in a cylindrical cell (diameter 
d = 44mm, length L = 60mm) and it is allowed to move 
along the optical axis. The turbid solution is made of 
an aqueous solution of poly-disperse silica particles (Lu- 
dox PW-50, average particle diameter ~ 50nm). Three 
volume fractions, </>i = 3.1 x 10 -3 , 02 = 7.8 x 10 -3 and 
03 = 15.6 x 10 -3 were used, with corresponding trans- 
port mean free paths Z* ~ 17.3mm, ZJ ~ 6.7mm and 
Z3 ~ 3.8mm. The light transmitted by the object and 
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FIG. 1. (Color online) Setup for backscattering DGI. The 
zoom shows a detail of the scattering cell with the six photo- 
diodes used for the bucket detection. 
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FIG. 2. (Color online) Setup for standard imaging. The ob- 
ject is back illuminated with a ring of speckled light shaped 
by using the lens Li, a diaphragm and a stopper. The lens 
Z/2 realizes a 1:1 imaging of the object on the CCD sensor. 

backscattered by the medium is collected by six photodi- 
odes placed in a ring configuration (ring radius = 15mm) 
outside the cell around the object. Such a configuration 
ensures that the average output signal from the six pho- 
todiodes can be used as an effective bucket detector. In- 
deed, the transport mean free paths of our medium are 
much smaller than the average contour length L c that 
photons travel from the injection point to the escaping 
point at the photodiodes positions. Thus, thanks to the 
backscattering detection scheme, the light reaching the 
photodiodes is completely randomized and the measured 
signal is proportional to the overall power injected into 
the solution and transmitted by the object. This implies 
that, in a blank measurement with no object, Si lk = a 62, 
where £2 = f A h{^2 )^x 2 and a is a factor which takes 
into account any unbalancing (beam splitter, detectors, 
random medium) between the two arms. 

The ghost imaging data analysis is carried out by using 
the DGI algorithm [7] based on the measurement of the 
observable 

<0_(x 2 )> = (5 1 / 2 (x 2 )) - |g(5 2 7 2 (x 2 )) (1) 

where 5i and 62 are the bucket signals collected in the 
object and reference arms, respectively, and ( • • • ) is per- 
formed over independent speckles configurations. From 
Eq. (pp) we measure the fluctuations of the transmission 
function 5T m (x), while the measured spatial average of 
the transmission function is T m = (Si)/(a(S2))- The 
measured transmittance of the object is thus computed 
asT m (x) =T m + JT m (x). 

The recovered ghost images are compared with stan- 
dard imaging measurements performed with the setup 
shown in Fig. [2] This setup is in a sense the reverse 
of the GI one since illumination is performed backward 
and detection forward, and is similar to the ones used 
in biomedical tissue imaging. A ring of speckled light, 
formed by reshaping the speckle beam with a diaphragm 
and a stopper, diffusively illuminates the object from the 
back satisfying the turbid medium condition I* <C L c . A 
macro objective (Nikon AF Micro Nikkor 60mm f/2.8D) 
realizes the imaging of the object onto the CCD sensor 




2 4 x (mm) 6 



FIG. 3. (Color online) Images of a thin black cardboard (sec- 
tion 1.8mm x 8mm, thickness s c± 400 /xm) inside a turbid 
solution with I* — 17.3mm at hi — 3.3mm (a), hi — 7.3mm 
(b) and /13 = 11.3mm (c), obtained with DGI (left column) 
and standard imaging (central column). The right column 
plots the corresponding horizontal sections averaged over the 
vertical dimension of the image: DGI (red curve), standard 
imaging (black curve), theoretical model (blue dashed line). 

with a 1:1 magnification. 

In our experiments we considered two simple objects 
characterized by a binary transmission function, with 
T(x) = 0, 1. The first object was a thin black cardboard 
of section 1.8mm x 8mm and thickness s ~ 400 /im. The 
turbid solution with I* — 17.3 mm was used. Fig. Ore- 
ports three examples of images recovered via DGI (left 
column) and standard imaging (central column), to- 
gether with their horizontal sections averaged over the 
vertical dimension of the image (right column). The 
DGI images were obtained by averaging 6000 indepen- 
dent speckle configurations. The figure shows that, as 
h is increased, the visibility of recovered images (both 
standard and DGI) becomes smaller because the central 
part of the image, where the object is totally absorb- 
ing (T = 0), becomes increasingly transmissive, passing 
from T rn ~ 0.1 (h = 3.3mm) to T rn ~ 0.5 (ft = 11.3mm). 
The figure shows also that the matching between DGI 
and standard imaging is excellent, although, as expected, 
DGI suffers of a much lower SNR. The latter one can be 
easily improved [7[ by increasing the number of measure- 
ments. The blue dashed lines in the right columns are 
the result of a simple model for DGI described below. 

Under the assumption of a thin object of thickness s, 
located at a depth ft from the entrance face of the cell 
such that s <C ft <C /*, the light travelling along the 
distance ft can be considered as undergoing only single 
scattering events. Thus, the light that hits the object 
(as sketched in Fig. 2(A)) is made of two main contribu- 



FIG. 4. (Color online) (A): backscattering scheme showing 
that the bucket signal Si is composed by three main contri- 
butions: (a) straight light hitting the object at a depth h, 

(b) forward scattered light that illuminates the object and 

(c) scattered light collected without intercepting the object. 
(B): scheme for the fraction sk of the light scattered in the 
forward direction that hits the object. 



tions: (|2a|) the straight non-scattered light that reaches 
the object with a probability 0h = exp(— h/l*) given by 
the Lambert-Beer (L-B) law Q and (|2EJ) the light that, 
after being scattered with probability 1 — fa, reaches the 
object with probability Sh determined by geometrical fac- 
tors. Hence, we may write this intensity as 



/?/Ji(xi) + 



+ (1 - Ph)e h 



A h 



=( x l) 



(2a) 
(2b) 



where J A 2i(x / 1 )dx / 1/ A4& is the incident average inten- 
sity and i sc ( x i) is the distribution of the scattered 
light, totally uncorrelated to 2i(xi), normalized so that 
J Ab (z ac (xi))dxi/Afe = 1. The bucket signal Si is given 

by J A 2-[^(xi)T(xi)<2xi plus a third contribution com- 
ing from the scattered light that does not pass through 
the object. Si can be written in the following way 



Si {fa I T(xi)/i(xi)dxi + 

S Ah h{A)dA 



(3a) 



{l ~ Ph)£h A b J Al 

(l-fa)(l-s h ) [ /i(xi)dxi) 

J Ah } 



(xi)T(xi)dxi+ (3b) 



(3c) 



Note that if the object is placed on the surface (h = 0), 
fa = 1, we recover the common definition of the bucket 
signal ([3a]) used in absence of the turbid medium, regard- 
less of eh- The dimensionless factor en reads 



£h 



hj 



dz fa w(z, h) 



(4) 



which is an average along the object depth h of the prob- 
ability u(z,h) that scattered light, at a position z from 
the surface of the cell, hits the object, weighted by the L- 
B factor fa (see Fig. HfB)). This probability corresponds 




FIG. 5. (Color online) Plots of the experimental data (sym- 
bols, whose dimensions are equivalent to error bars) and the- 
ory (solid curves) for the absorbing region of the object T m as 
a function of its depth h for three different turbid solutions. 



to the fraction of light scattered within a maximum solid 
angle Q(z) subtended by the object (with area Ab), nor- 
malized to ^ = J 4?r sin 2 cj)dQ. We assume Rayleigh scat- 
tering with an incident polarized electric field that forms 
an angle <fi with the scattering direction. The integral in 
Eq. (J2J) is computed numerically. 

Combining Eq. (pQ) and Eqs. (j3j), with the assump- 
tion of uniform illumination ((ii(xi)) = (2i), (22 (X2)) = 
(22)), and taking into account that 22 (x 2 ) and i 3C (x.i) 
are uncorrelated ((2 2 (x 2 )i sc (xi)) = (2 2 (x 2 ))(i sc (xi)) = 
(2 2 )), we derive an expression for the measured trans- 
mittance of the object T m (x) in terms of the real T(x), 
which reads 

T m (x) = ^T(x) + (1 - fa) \e h T + l-e h ] (5) 

where T = J A ^ T(x.i)dx.±/Ai J is the spatially averaged 
transmittance of the object. As expected, Eq. (j5j) pre- 
dicts that, in the case of non-turbid media or in the 
case of objects placed at the surface of the cell, when- 
ever h/l* ->> 0, T m (x) -» T(x). But remarkably, al- 
though based on the assumption that h <C Z*, Eq. (|5j) 
predicts also the correct behavior of T m (x) for highly 
turbid media or objects deeply inside the scattering cell 
(h/l* — >> 00), for which fa —> and Sh 0. In these 
cases, indeed, the object becomes invisible and, consis- 
tently, Eq. (j5j) predicts T m (x) -> 1. When Eq. (j5j) is ap- 
plied to analysis of the images of Fig. [3] (blue dashed lines 
in the third columns), the agreement with the experimen- 
tal data is excellent in correspondence of the absorbing 
zones (T(x) = 0) of the object, while is somewhat less 
accurate for the transmissive zones (T(x) = 1). Over- 
all, the simple model of Eqs. (j2|3p is able to capture the 
essential physics of DGI in turbid media. A more quanti- 
tative analysis of the data of Fig. [3] is reported in Fig. [5j 
where we compare, as a function of h, the behaviors of 
the expected (Eq. flfaj)) values T m of the absorbing zone 
(T(x) = 0, solid curves) with the experimental data, for 
the three solutions with I* ~ 17.3mm, /J 6.7mm and 
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FIG. 6. (Color online) Images of a small absorbing sphere 
(diam = 900 /mi) inside a turbid solution with I* = 6.7mm at 
hi — 2.4mm (a), hi — 4.4mm (b) and h% — 6.4mm (c), ob- 
tained with DGI (left column) and standard imaging (central 
column). The right column plots the radial profiles averaged 
over the azimuthal angle: DGI (red curve), standard imaging 
(black curve), theoretical model (blue dashed line). 



in 



3.8mm. The agreement between theory and ex- 



periment is quite good for the I* = 17.3mm curve but 
becomes less accurate at higher turbidities, where the 
presence of increased multiple scattering reduces the va- 
lidity of the assumptions used in the model of Eqs. (j2!3p . 

Our results were also validated by measuring an ab- 
sorbing sphere much smaller (diam = 0.9mm) than the 
beam area (T ~ 0.99) in a turbid solution with P = 
6.7mm. Figure [6] reports the images retrieved with DGI 
(left column) and standard imaging (central column), to- 
gether with their corresponding radial profiles (right col- 
umn) obtained averaging the images over the azimuthal 
angle. We notice that, as for Fig. [3j the object becomes 
less visible as h is increased. In correspondence to the 
absorbing region of the object, we obtain values that 
ranges from T rn ~ 0.30 (hi = 2.4mm) to T rn ~ 0.61 
(hs = 6.4mm), a result that is equivalent for both DGI 
and standard imaging. The agreement between experi- 
mental results and the theoretical model is excellent, as 
shown in the third column of the figure. 

In this letter we have shown that DGI can be profitably 
used in a backscattering configuration for the imaging of 
small absorbing objects immersed in a turbid medium, 
in proximity of its surface. Spurred by the recent debate 
about the potentiality of GI for the imaging of objects 
in the presence of turbulence or scattering [11], we have 
quantitatively compared DGI with a standard imaging 
method. Our results show that the two techniques per- 
form almost identically and are equally affected by the 
presence of multiple scattering when the object is deeply 



immersed in the medium (h > I*). This feature demon- 
strates that GI is not immune from multiple scattering, 
exactly as it happens for GI when there is turbulence 
between the beam splitter and the object or the CCD 
[Tlj . However, there are situations where backscatter- 
ing DGI may turn out to be very convenient, such as 
for example in biomedical tissue imaging for the early 
detection of pigmented skin lesions. In these cases, the 
existing optical methods [lij are either too qualitative 
(such as epiluminescence imaging or dermoscopy [16|) or 
rather complex and expensive as the ones based on Opti- 
cal Coherence- [TtJ and Diffuse Optical- Tomography |18| 
or Diffuse Reflection-correlation spectroscopy [19]. We 
therefore believe that backscattering DGI has the poten- 
tialities to become, in the next future, a valid imaging 
tool alternative or complementary to the current state of 
the art imaging techniques. 
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